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ABSTRACT 

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show 
slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 
to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may 
result from scattered and thermal light echoes off of extended circumstellar material. We present a 
semi-analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to 
provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous 
dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar 
shell 14-15 lt-yr from the progenitor, and containing < 0.02 Mq of carbon-rich dust, can explain 
the broadband spectral and temporal evolution. The size, mass and dust composition are in good 
agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin 
of slowly-changing high- velocity [O I] and Ha lines is also considered. We propose an origin in shocked 
high- velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving 
circumstellar material, as is the case with hot spots in SN 1987A. 

Subject headings: supernovae: individual (SN 1980K) — dust, extinction — circumstellar matter — 
methods: numerical — radiative transfer 



1. INTRODUCTION 

Nine supernovae (SNe) have been observed in the spi- 
ral ga laxy NGC 6946 (d = 5.9 Mpc, Karachcntsc v et all 
2000) in the last century, earning this nearby super- 
nova factory the nickname "The Fireworks Galaxy" and 
making it a virtual laboratory for studying a variety of 
SN characteristics from progenitor identification through 
evolution into early remnant stages. With the Hubble 
and Spitzer Space Telescopes (HST and Spitzer, respec- 
tively) and large ground-based observatories all available 
during the last decade, a great deal of understandable 
interest has been paid to its recent SNe 2002hh, 2004et, 
and 2008S. Continued monitoring of its older population 
(SNe 1917A, 1939C, 1948B, 1968D, 1969P, and 1980K) 
provide an extraordinary opportunity to connect the 
early evolution (< 10 yrs) of young SNe with evolved 
remnants (e.g. Cas A, N 103B). In this paper, we re- 
port and analyze optical and infrared (IR) data collected 
between 2004 and 2010 of SN 1980K. 
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First discovered on 1980 Oct 28 (|Wild fc Barbonl 
1980), this core-collapse, linear (or type II-L) SN was re- 
covered in Ha imaging in 1987 and spectroscopy in 1988 
(|Fesen fe B ecker 1990), which revealed a faint continuum 
underlying a complex of high- velocity Ha and [O I] lines 
that had not been present during the fi rst two years 
(|Barbieri et all Il982t jUomoto fe Kirshnerl 11986ft . Con- 
tinue d observations ( Uomotol 119911 : IFesen fc Matonickl 
119941: IFesen et all I1995L 119991 ) have shown the total 
fluxes to have declined by about 25% over more than 
a decade, with the centers and fluxes of the individual 
high- velocity components changing on timescales of a few 
yea rs. A detail e d con sideration of possible mechanisms 
led IFesen et al.l (|1999( ) to conclude that these lines orig- 
inate in shocks driven by the SN blast into pre-existing 
circumstellar material (CSM). The presence of such ma- 
terial has been inferred both f rom IR excesses o bserved 
up to o ne ye ar after outburst ([Dwek et al.lll983l ). which 
iDwekl (|1983l ) argued arose from thermal echoes off of 
nearby CSM, and from 6 and 20 cm radio evolution 
(jWeiler et al.lll992l ) consistent with the forward blast im- 
pacting dense CSM laid down by the progenitor red su- 
pergiant a few tens of thousands of years prior to core 
collapse. 

The presence of a faint continuum within all late-time 
spectra, and the slow monotonic decrease in flux in all op- 
tical waveband s, can also arise from a scattered-light echo 
(|Couderdll9"39l ) illuminating a large CSM shell or inter- 
stellar dust sheet. Whether from scattered or thermally- 
reprocessed SN light, an echo arises when dust redirects 
SN light back into the line of sight, which arives at the 
observer later than the original signal by the light-travel 
time along the longer path length. Indeed, a scattered 
echo h as been considered for SN 1980K by lChevalie"rl 
(1986), who concluded that the flux evolution through 
the first ~ 1 year was consistent with a scattering-echo 
model, but more observations were necessary to distin- 
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Figure 1. Tiling of images of SN 1980K from optical through far-infrared. Note that the image scale (arcsec) changes between the top 
and bottom two rows, (a) HST/WFPC2 F606W on 19 Jan 2008. (6) Gemini r' on 05 Aug 2005. (c) Spitzer/IRAC 4.5 fim on 05 Jan 2011. 
(d) Spitzer/IRAC 8.0 /im on 18 Jul 2008. (e) Spitzer/MTPS 24 fim on 29 Jul 2008. (f-i) Herschel/PACS and SPIRE on 03 Oct 2010 at 70 
fim, 100 fim, 160 fim, and 250 fim, respectively. 



guish it from radioactive-energy deposition. 

In this paper, we present optical photometry of SN 
1980K made by the HST and Gemini North, and mid-IR 
photometry made by Spitzer, all between 2004-2010, as 
well as upper limits in the far-IR from Herschel Space 
Observatory imaging made in 2010 (©. Combined with 
previously-reported optical photometry, we show that 
the last decade of evolution is consistent with a slow, 
monotonic fading across all wavebands ($3]) that can be 
explained by scattered and thermal echoes from a large 
(r ~ 14 lt-yr) shell of circumstellar dust. The basic be- 
havior and modeling of light echoes are discussed in <2J 
followed by an exploration of the parameter space of pos- 
sible shell properties, as well as a new, semi-analytic dust 
radiative-transfer model, in SjXJ The consistency of the 
proposed light-echoes with previous analyses of the high- 
velocity emission-line structures is discussed in £j5j 

2. OBSERVATIONS 
2.1. Imaging 

In the mid-infrared, SN 1980K has been imaged at least 
yearly between 2004 Jun and 2008 Jul with Spitzer in 
the same field as SN 2004et with both the InfraRed Ar- 
ray Camera (IRAC) and Multiband Imaging Photome- 



ter for Spitzer (MIPS) as part of the original cold mis- 
sion (GO-159, 20256, 20320, 30494, 40010), and con- 
tinuing through 2010 with IRAC as part of its ongo- 
ing warm mission (GO-60071, 70008), often as part of 
the ongoing work of the SEEDS collaboration (Search 
for Evolution of Emission from Dust in Supernovae, 
P.I. M. Barlow). Pipeline-calibrated imaging data were 
retriev ed from the Spitzer arc hive and processed with 
mopex (jMakovoz k, Khan 2 0"05() to achieve enhanced res- 
olutions of 0'.'75 pix" 1 for IRAC and 1"5 pix" 1 for MIPS. 

SN 1980K was imaged once on chip 2 of the HST 
Wide Field and Planetary Camera 2 (WFPC2) in 
the F606W and F814W filters on 2008 Jan 19 (GO- 
11229). Pipel ine-reduced data were combined using 
multidrizzle (Kockc moer et aLll2002[ ) to achive a final 
pixel scale of O'.'l pix" 1 . Two epochs of optical observa- 
tions were also made with Gemini North in g'r'i filters 
on 2005 Aug 05 and 2006 Jul 17 using the Gemini Multi- 
Object Spectrograph (GMOS) in imaging mode. Images 
were reduced and stacked using the IRAF gemini pack- 
age, and instrumental counts were calibrated to Johnson- 
Cousins VR r I r ma gnitudes using transformations from 
IWelch et all (l2007n . 

NGC 6946 was imaged in the far-infrared by Her- 
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schel in all wavebands on 2010 Mar 10 as part of the 
KINGFISH program (R. Kennicutt, P.I.). One-arcmin 
fields surrounding the position of SN 1980K were ex- 
tracted from the enhanced Photoconductor Array Cam- 
era and Spectrometer (PACS) and Spectral and Photo- 
metric Im aging REceiver (SPIRE ) data products as ex- 
plained in Kennicutt et al. (2011). Images of the field at 
representative wavelengths are shown in Fig. [TJ 

2.2. Archival Spectra 

Optical spectra of various wavelength ranges and res- 
olution have been published for SN 1980 K from max- 
lmum light in 1980 through late 1997 (iBarbon et al l 



1 1982. lUomotq fe Kirshnei1ll986|: ILeibundeut et al.lll991l 
Uomotol 11991. IFesen fc Beckerl Il990l : IFesen fc Matonickl 
Il994t IFesen et al.l Tl995. 199j|). Fo r this work, optica l 



1982) 



1991) 



spectra from the first ~ 100 days (IBarbon et al.l 
that were digitized and flux calibrated by iBenettil 
were acqu ired from the Asiago Supern ov a Archive. Spec- 
tra from lUomoto fc Kirshnerl J 19860 . IFesen fc Beckerl 
(fl990T) . and IFesen et all (11999) were extracted using 
DEXTEfF 1 ] from fully-printable scans or electronic articles 
archived in the Astrophysical Da ta Syste nP^I. Th e late- 
time spectrum from 1988 Aug (jFesen fc Beckerl I1990D 
was integrated over the Johnson V and Cousins R c fil- 
ter response funct ions, and the spectrum from 1997 Nov 
(jFesen et al. 1999) over V, R c and I c , yielding estimates 
of V = 22.9 and R c = 22.1 in 1988, and V = 23.3, R c = 
22.4 and I c < 22.9 in 1997. These are consistent with the 
broadband magnitudes V = 2 2.8 ± 0.2, R = 21. 9 ± 0.1 
and / = 22.2 ± 0.3 reported by ILeibundeut et al.l (|1993f ) 
between 1990-1992. 

3. PHOTOMETRIC ANALYSIS 

3.1. Astrometry 

Given that SN 1980K is >25 years past maximum light 
in our data, it is important to verify the position of 
the point source identified as the SN. 15 stars in the 
HST/WFPC2 field were extracted from the USNO-B1 
astrometric catalog (| Monet et al.l I2003D and compared 
to positions in the HST and Spitzer data to calibrate the 
world-coordinate systems using standard IRAF routines. 
These show that the point source at the center of panels 
(a-e) in Fig. [J are within O'.'l of the VLA position of the 
SN as given in iMontes et al.l (|1998| ). In particular, SN 
1980K is the right-hand point source at the very center 
of the HST image shown in Fig. [Da. 

3.2. Photometry 

In optical and mid-IR imaging, photometry was per- 
formed using the daophot PSF-fitting routines within 
IRAF, which also provide a robust estimate of uncer- 
tainty using a combination of Poisson noise, flat-fielding 
and absolute calibration errors. The resulting fluxes and 
errors, as well as earlier data estimated from archival 
spectra, are shown in Fig. __] Also shown in grey are 
ranges of slopes that fit the data with standard least- 
squares. The general trend is for a slow fading in time in 
most all wavebands, although most data are consistent 
with the flux remaining approximately constant over the 

11 http:/ /dc. sah.uni-heidelberg.de/sdcxter 

12 http://adsabs.harvard.edu 



time period sampled. For this reason, the weighted aver- 
age values for each filter are presented in Table Q] along 
with the least-squares slopes, and the averaged spectral- 
energy distribution (SED) is plotted in Fig. H 

No point sources were detected at the position of the 
SN in Herschel data. Although there is a point-source 
like feature at the correct position in the 160 /im frame 
(Fig. [lji) , it has the same brightness as horizontal and 
vertical banding throughout the image, suggesting it is 
coincident noise. To estimate the upper limits in each 
frame, a grid of PSFs were added to each image with 
varying brightness until daophot could detect them. The 
resulting PSF fluxes are close to the upper limits qu oted 
for each imager by the Herschel Science Central, and 
have been listed in Table [T] and plotted as downward 
arrows in Fig. [3] 

3.3. Interpretation 

The mid-IR SED shown in Fig. [3] is consistent with 
thermal emission from warm dust. Regardless of its 
origin, a zeroth-ordcr estimate of the composition and 
mass of emitting particles can be made by assuming that 
all grains are visible (i.e. negligible internal extinction, 
which is reasonable for mid-IR wavelengths) , have a sin- 
gle temperature T, and emit as blackbodies modified 
by their mass-absorption coefficients k\ — ira 2 Q a b s /Tn g 
where a is the grain radius, Q a bs is the absorption effi- 
ciency and m g is the grain mass. The emerging spectrum 
is then given by 



Fx(T) = 



kxB x (T)M 
D 2 



(1) 



where D is the distance to the dust and M is the 
total dust mass (|Dotv fc Leund U99l . Three com- 
positions were tested using the grain propertie s from 
iWeingartner fc Draind (|2001| ) and lLaor fc Draind (j!993l ) 
- pure astronomical silicate (Si), pure graphitic solid (C), 
and an even mixture of both (C+Si) - with a standard 
a -3.5 distribution of grain sizes (iMathis et alj|T977l here- 
after MRN). Note that at mid-IR wavelengths, the ab- 
sorption coefficient n is relatively insensitive to grain size. 
A large range of dust temperatures were sampled for each 
composition, with the dust mass at a given temperature 
chosen to minimize the x 2 residual between the photo- 
metric data and synthetic fluxes found by integrating the 
model SED over the Spitzer and Herschel filter profiles. 

The best-fit models are shown in Fig. (3]a, representing 
2 x 10~ 4 Mq of Si dust at 200 K, 8 x 10" 4 M of C dust at 
230 K, and 2 x 10~ 4 M of C+Si grains at 245 K. While 
the C and C+Si models come close to matching the mid- 
IR SED, the fits differ from the data by x 2 > 30 suggest- 
ing that a single-temperature model is insufficient. We 
repeated the same exercise using two dust components, 
where the best-fit mass for each model was determined 
using the downhill-simplex algorithm amoeba( Press et al.l 
1992). For a given temperature of the warm-component, 
a large range of cool-dust masses can fit the SED, due 
to the high upper limits at far-IR wavelengths. Fig. __}? 
shows the best-fit models using the lowest masses for each 
cool-dust component. While the C and C+Si models 
have x 2 ~ 1, the Si models all have x 2 > 21, suggesting 

13 http:/ /herschel. esac.esa.int/science_instruments.shtml 
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Table 1 

Photometric Behavior of SN 1980K 



Ao 


Avg. Flux 


Best-Fit Slope 


N 1 




Q 


F(\) at 23 yr 


F(X) at 28 yr 


w 


(pJy) 


(/iJy vr" 1 ) 






0*Jy) 


OJy) 


0.55 


1.7±0.4 


-0.044 ± 0.039 


5 


0.873 


1.8±0.4 


1.6±0.4 


0.64 


3.1 ±0.7 


-0.043 ± 0.057 


1 


0.091 


3.0±0.7 


2.8±0.7 


0.80 


1.9 ±0.5 


-0.006 ± 0.100 


1 


0.740 


1.9±0.7 


1.9±0.7 


3.6 


1.8 ±0.5 


-0.067 ± 0.103 


9 


0.831 


2.1±0.7 


1.8±0.7 


4.5 


5.9 ±0.6 


-0.352 ± 0.090 


12 


0.024 


7.7±0.7 


5.9±0.8 


5.8 


16.8 ±5.2 


-1.936 ± 1.769 


5 


0.968 


22±10 


13±10 


8.0 


51.9 ±13.4 


0.706 ± 3.147 


11 


0.552 


50±21 


53±21 


24. 


460. ± 28. 


-14.73 ± 8.15 


10 


0.138 


503±50 


430±50 


70. 


< 5900. 




1 




< 5900. 


< 5900. 


100. 


< 5900. 




1 




< 5900. 


< 5900. 


160. 


< 9900. 




1 




< 9900. 


< 9900. 


250. 


< 15000. 




1 




< 15000. 


< 15000. 



1 Number of epochs in fit and average 

2 The incomplete gamma function Q gives the goodness-of-fit probability. 



that the dust responsible for the mid-IR emission must 
have some carbonaceous component. 

The dust masses presented in Table [2] represent the 
instantaneous amount of dust necessary to produce the 
average SED. Since dust cooling times are short com- 
pared to the slow fading of the light curves, we expect 
that either a small amount of dust is being heated by 
a source with decreasing energy, or the volume of emit- 
ting dust is changing with time. In the former case, the 
dust mass present in the system is roughly that found 
in this section, while in the latter, the total dust mass 
is much larger than that which is emitting at any given 
time. Light echoes are one natural mechanism for illu- 
minating a slowly-changing volume of dust with time, 
where for our timescales involved, the initial SN light 
pulse would be heating and scattering off of dust in a 
large, pre-existing circumstellar shell. In the next sec- 
tion, we present a new method for modeling such echoes, 
and explore whether such echoes can reproduce the late- 
time data. 

4. LIGHT ECHOES 

4.1. Formalism 

Light echoes occur when a pulse of light from a variable 
object (e.g. a SN) interacts with dust and arrives at Earth 
some time t after the pulse is directly observed. Optical 
light that scatters from that dust into the line of sight 
will be observed as a scattered-light echo, while prompt 
thermal re-emission from dust heated by the pulse will 
be observed as an infrared (or thermal) echo. Scattered 
echoes from SNe h ave b een d iscusse d in e .g. lChevalie"r1 
(p8l) . ISugermanl (|2003l ). and iPatatl (120051) whi le ther- 
mal echoes h ave be en presented in|Bode fc Evans! (|1980t ). 
[Dwell (|l98l . and lEmmering fc Chevalierl (|1989ft . The 
following is a brief synopsis of the relevant properties of 
light echoes. 

An echo lies on the locus of points equidistant in light- 
travel from the source and observer, i.e. an ellipsoid with 
the source and observer at its foci. In the neighborhood 
of the source, this ellipsoid can be approximated by the 
"echo equation" 

2 P 2 & 
Z ~ 2ct ~2' 



(2) 



Dust being illuminated at a time t after the SN is directly 
observed, and located a distance r from the source at an 



angle 9 from the line of sight will have line-of-sight depth 
z = r cos 9 and projected distance p = r sin 9 on the plane 
of the sky. We adopt distances and times in lt-yr and yr, 
respectively, which simplifies this eq uation since c = 1 . 

Proceeding as in ISugermanl (|2003t l. the flux scattered 
off one dust grain of radius a at position r from the SN 

dF sca (X, r, a) = — 2 (3) 

where C — -ko?Q and Q sca is the grain scattering effi- 
ciency, F(X) is the fluence (i.e. time-integrated flux) at 
the surface of the SN, and <E> is the scattering phase func- 
tion for a given scatte r ing an gle /i = cos 9. We adopt the 
iHenvev &: Greenst ein (194l|) phase function 



$(^,A,a) 



l-.9 2 (A,a) 



[I + .g 2 (A,a)-2g(A,a) M ] 3/2 



(4) 



with <7(A, a) measuring the degree of forward scattering 
for a given grain. The total flux F sca integrated over 
the duration At of the outburst from a single scatter- 
ing is found by multiplying Eq. (j3|) by the dust density 
nd(r,a), integrating over the scattering volume and all 
grain sizes. If the optical depth of the cloud is low, then 
one can adopt the single-scatte ring plus attenuation ap- 
proximation (SSA, Patat 2005), in which the flux that 
reaches the dust grain is extinguished by e~ Tr along the 
radial path from the source to the grain, and further by 
e~ Tz out of the dusty medium along the line of sight 
to the observer. Implicit in this approximation is that 
no photons are scattered into the line of sight by other 
grains of dust within the medium. If the medium con- 
sists of a distribution /(a) of different grain sizes such 
that n c i(r,a) = n(r)f(a), the total scattered flux from a 
volume element d 3 r arriving at Earth will be 



F M 



1 (A,r) = 

C sca (A,a)F(A)$(^,A,a)n /( a ) 
l§Ti 2 D 2 r 2 



(5) 



^— ry— t 



dad 3 r 



where n(r) — no(ro/r) p allows for a power-law depen- 
dence of density with radius, and D is the distance 
to the SN. Examples of grain-size distributions include 
the previously-mentioned M RN function and those of 
iWeingartner fc Drainej (2001). 
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Table 2 

Results of two-component dust modeling 



Warm Component Hottest Cool dust Coldest Cool dust 

Model T warm (K) M (M ) T cool (K) M (M ) T cool (K) M (Mq) 

C 300 - 330 4 x 10" 5 180 10" 3 55 3. 

Si 450 - 480 10" 4 140 2 x 10" 3 45 1. 

C+Si 330-350 10~ 5 160 lO" 3 50 3. 



A grain of radius a heated to temperature T will radi- 

F(X) = C ahs (X,a)B x (T) (6) 

where C a b s is the grain absorption (and emission) cross 
section and B\ (T) is the Planck blackbody function. By 
analogy with Eq. ([5]), the total thermal emission from all 
grains within a volume element arriving at earth will be 

F th (X,r) = J — 3 (-j e ~dadr 

(7) 

4.2. Illuminating Spectra 

The optical spectrum F(X) that produces scattered- 
light echoes is taken to be the time-integrated fluence 
over the first ~ 100 days, constructed by interpolat- 
ing each spectrum onto a regular wavelength grid and 
integrating them over time using the trapezoidal rule. 
When considering only scattered-light, one may model 
the echoes using the observed (i.e. reddened) spectrum 
since it has suffered the same interstellar extinction as 
any echo flux, however as we wish to also model the ther- 
mal response to the illu minating fluence, th e dereddened 
spectrum must be used. iFabbri e t al. (2011) adopt a red- 
dening of E(B-V) = 0.41 toward the nearby SN 2004et, 
noting that the foreground Galactic contribution is esti- 
mated to be E(B - V) = 0.34. If SN 1980K is located 
deep within or behind the disk of NGC 6946, one may 
expect it to have produced larger-radius light echoes, 
as have been o bserved f or example around SNe 1 987A 
(IXu et al.lll995h. 1993.T (ISugerman fc Crott3l200l . and 
2003gd (|Sugermanll2005l ). Careful inspection of all opti- 
cal images show no such echo features, thus in our mod- 
eling, we consider that while SN 1980K may suffer the 
same extinction as SN 2004et, it may also lie on the near- 
side of its host galaxy and its extinction may be closer 
to the lower, Galactic value. In what follows, we desig- 
nate these the "high" and "low" extinction models, re- 
spectively. The fluence spectrum, dereddened by the full 
E(B - V) = 0.41 is shown in Fig. gja. Note that since 
this fluence spectrum ends at 7500 A, it is not possible 
to model the J-band scattered flux of any echoes. 

Modeling the thermal echo of a SN requires considering 
all the photon heating sources of the dust, which includes 
not only the long-duration optical pulse but the preced- 
ing short-duration but highly-ener getic UV puls e, whic h 
has been variously m o deled by iGezari et al.l $2 008), 
Chevalier fc Franssonl (12008ft . iNakar fc Saril (|2010t l. and 
Rabinak fc Waxmanl (|2011| ). The blackbody parameters 
presented in each paper were used to create fluence spec- 
tra over th e first ~ 2.3 days, as shown in Fig. Al- 
though the Cheva lier fc Franssonl (|2008f) model was pro- 
posed for Ib/c SNe, its fluence spectrum describes a rea- 
sonable average between the other three models consid- 



ered, therefore we adopt this spectrum as describing the 
initial UV pulse. 

4.3. Constraints 

Fig. [TJa shows that there is no resolved flux around 
SN 1980K, and examination of the central source shows 
that it is consistent with a point source, thus any echo 
must come from material such that its flux remains un- 
resolved. Experimentation with marginally-resolved ar- 
tificial sources show that an echo can remain unresolved 
with a projected distance p up to 1.2 times the full-width 
at half-maximum (FWHM) of the PSF, which for the 
WF2 chip yields p < 0.2'.' or roughly 19 lt-yr. At t = 28 
yrs, an echo would be at z ~ — 7 lt-yr and r ~ 20 lt-yr. 
For an echo to be present t yrs after outburst, material 
must be at least t/2 lt-yr from the source. Since there is 
no other spatial information available, we will model the 
echo as arising from a simple sp herical shell with outer 
radius between 14 and 20 lt-yr. Chevali er fc EmmerineJ 
(1989) argue that red supergiant winds expanding into 
ambient interstellar media will have radius of order 16 
lt-yr, consistent with these constraints. 

Since the average light curves presented in Fig. [2] arc 
monotonically fading in almost all wavebands, we opt to 
model the system using the average optical and mid-IR 
data from 23 and 28 yrs after outburst, which cover the 
period of multi-waveband Spitzer observations. Fluxes 
and their associated uncertainties for each epoch were 
taken as the least-squares fits to each monochromatic 
light curve, and are listed in the two rightmost columns 
in Table [TJ Lacking any temporal data in the far-IR, 
we assume that the upper limits from Herschel apply 
throughout the epochs studied. 

4.4. Modeling Echoes 

The optical and thermal echo SED resulting from a 
SN will be the sum of the light scattered from the short 
(~ 2 day) UV and long (~ 150 day) optical pulses, as 
well as the thermal emission resulting from the dust be- 
ing heated by those two input spectra. The UV pulse 
contributes negligibly to the optical echo, however both 
spectra are important sources of heating for thermal 
emission. 

The equilibrium temperature of dust illuminated by 
spectrum F(X) (and excluding collisional heating, which 
is reasonable in low-density, slow-m oving, unshocked 
circu mstellar media) is given by (e.g. iRowan-Robinsonl 
fl980h 

4 J C ahs {X,a)B x (T)dX=^ J F(X)C ahs (X,a)e- T < dX. 

(8) 

Numerically, the dust temperature for a given grain size 
a can be calculated using a root-finding algorithm such 
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Figure 2. Optical (top) and mid-IR (bottom) light curves of SN 
1980K. Wavebands are indicated in the top right of each panel. 
Grey regions indicate the range of slopes consistent with the data 
using least-squares fitting. 
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Figure 3. Average SED of SN 1980K using data in TableQ] also 
showing model spectra for simple dust emissio n mo dels using (a) 
one and (b) two temperature components (see £|3,3I1 . Dotted lines 
denote silicate dust, dashed lines denote graphitic dust, and solid 
lines show an equal mixture of both. 



as Brent's method (jPress et al.lll992h which varies T un- 
til both sides of Eq. (JTJ) are equal. If dust is not heated 
beyond < 1000 K, the majority of its emission will be 
in the mid-IR and will contribute insignificantly to re- 
heating of neighboring dust. A thermal echo spectrum 
can thus be computed as the emission coming from the 
adjacent parabolic regions illuminated by the UV and op- 
tical pulses at a given time t after outburst by solving Eq. 
|8]) at each point along the echo parabola and summing 
the contributions using Eq. Similarly, the scattered 
optical component is computed by summing Eq. ([3]) in 
the same regions. 

The analytic expressions presented above are only ap- 
propriate in the optically-thin regime, since multiple 
scatterings within an optically-thick mediu m return pho - 
tons to a given path as well as remove them. Patat (2005) 
has shown that for a variety of dust geometries, the SSA 
breaks down for optical depths r > 0.1 and that a proper 
treatment requires Monte Carlo simulations. A variety of 
dust radiative-transfer Monte-Carlo models exist, how- 
ever they are numerically expensive, with the result that 
exploration of a large parameter space of dust geome- 
tries and densities to fit a given dataset can require tens 
of thousands of computing hours. A further complica- 
tion is that of resolution: most three-dimensional Monte- 
Carlo models track photon paths through a spatial grid 
of cells, making the run time and memory allocation de- 
pendent on N 3 where N is the number of cells in a given 
dimension. The UV pulse lasts only ~ 2 days, however a 
star's CSM can span tens of lt-yr. Without an adaptive 
mesh-refinement scheme, the spatial sampling needed to 
describe the dust's behavior in this narrow echo requires 
a prohibitively high numbers of grid cells. 
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Figure 4. SN Fluences used for modeling echoes from SN 1980K. (a) The composite optical spectrum of SN 1980K integrated over the 
first ~100 days, and deredde ned by E(B — V) = 0.41. (b) Different models for the shor t, prompt UV out burst. The solid black line is from 
Chevalier & Fransson (2008), while in grey, the dashed line is from Nakar & Sari (2010), dot-dashed from Rabinak & Waxman (2011), and 
dotted from Gczari ct al. (2008). 




Figure 5. Results of QuickSAND modeling of optical and thermal echoes to reproduce the average photometry of SN 1980K at 23 and 
28 yrs after outburst (Table [TJ. Each panel shows in greyscale the dust mass that best fits the data for a constant density (p = 0), 
low-extinction model with the given values of inner-shell radius ro, shell thickness b and carbonaceous-dust fraction c. The wedge at right 
shows the scaling between shading and the shell's dust-mass. Overplotted on each panel are contours of constant reduced x 2 at intervals 
of Xred = 0-5, with the outermost tracing Xred = ^- The results for high extinction are nearly indistinguishable, with the contour interval 
changed to unity and the outermost contour tracing Xred = ^- I n the neighborhood of the good fits, the p = 2 models are also nearly 
indistinguishable. 



It is for these reasons that we have developed the 
QuickSAND semi-analytical model presented in Ap- 
pendix[A] which approximates the results of Monte-Carlo 
dust radiative-transfer models in homogenous media to 
a high degree of accurac y. For an alternative semi- 
analytic model, please see iVarosi fc Dwe k ( 1999, here- 
after VD99). 

To test whether the proposed scattered and thermal 
echoes can reproduce the broadband photometry of SN 
1980K, we ran the QuickSAND model on a large parame- 
ter space of homogenous spherical shells, where each pos- 
sible shell is described by an inner radius ro, thickness 
6, density power-law index p f ^4. 1|) . and carbonaceous- 
dust fraction c (such that the amount of silicate dust is 
1 — c). Values of ro and b were tested every 0.5 lt-yr 
such that rO + b < 20 lt-yr (S3]), with p = and 2, 
and c — i, |, and 1 (recall modeling from £ )3. 31 showing 
silicate-dom inated dust does not reproduce the SEDs). 
As noted in iSugermanl (2003), one can use the observed 
spectrum when modeling a scattered-light echo since the 
observed spectrum has already suffered the appropriate 
extinction along the line-of-sight to the observer. How- 
ever, the extinguished spectrum is inappropriate for heat- 
ing of dust as the dust receives the unextinguished light. 
The QuickSAND SEDs were constructed using the un- 
extinguished optical and UV bursts in Fig. 0] and then 



reddened by E(B — V) = 0.42 to compare against data. 

Adopting the same grains as were used in £13. 3[ we used 
Brent's method to determine the dust density that mini- 
mizes the x 2 residuals between the QuickSAND model 
SEDs and data at both 23 and 28 yrs (Table [Q) for 
each set of input parameters. This is another advantage 
of our semi-analytic model over full radiative-transfer 
codes, since QuickSAND can automatically find the best- 
fit density rather than the user having to modify it by 
hand through trial and error. Another option is to al- 
low such models to find the best fit of all parameters, by 
using a multi-dimensional minimization scheme such as 
the downhill simplex, however there are many potential 
problems with this approach, namely, the minimization 
routine can settle into local minima, the output does not 
provide uncertanties or ranges of parameter space that 
fit data, and the runtime is significantly longer than step- 
ping through a grid of parameters while minimizing only 
one (e.g. density). 

The results for the low-extinction model are summa- 
rized graphically in Fig. [5j which shows in greyscale the 
dust mass that best fits the data for differing values of 
ro, &, and dust composition, assuming constant density 
with radius [p = 0). Overplotted in white are contours 
of reduced x 2 (i- e - X 2 P er degree of freedom) at inter- 
vals of Xred = 0.5 with the outermost countour tracing 
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Figure 6. Best-fit high (solid) and low-extinction (dotted) Quick- 
SAND models of the average optical and mid-IR fluxes at 23 (top) 
and 28 (bottom) years after outburst. Both models are for a shell 
with r =14 lt-yr, 6=1 lt-yr and c = 1.0, with M d = 1.5 X 10" 2 M 
and Xr e d = 1-6 m the high-extinction case, and = 1.8 X 10~ 2 
Mq and Xred = ^ or ^ ow extinction. 

Xred = ^- ^he best-fit models are thin shells with ro ~ 14 
lt-yr, 6=1 — 4 lt-yr, a medium dominated by carbona- 
ceous dust, an d a total dust mass of . 02 — .05 Mq. As 
normalized by Weinga rtner fe Drainel ()2001l) . the MRN 
grain distribution corresponds to a gas-to-dust mass ra- 
tio ~ 150, which yields a total shell mass ~ 3 Mq. 

Results for the higher-extinction model are nearly in- 
distinguishable, except the contours are spaced at inter- 
vals of Xred = 1- w hh the outermost contour tracing 
Xred = 4- F° r both low and high extinction cases, the 
p = 2 models also yield nearly indistinguishable results in 
the neighborhood of the good fits. In all models, the V- 
band optical depths of the good fits are < 10 -3 , i.e. the 
proposed CSM is optically thin. The SEDs for the best- 
fit low and high-extintction models are plotted against 
the data in Fig. [5] 

Such a structure may be the contact discontinuity be- 
tween the star's winds and the ambient ISM (also some- 
times called a wind-blown bubble). The size scale, mass 
and dust-grain properties are roughly consistent with the 
contact discontinuity a round SN 1987A reconstructed by 
iSueerman et al.l (|2Q05[ ). which was found to be a pinched 
prolate spheroid, with polar and equatorial axes of 20 and 
11 lt-yr, shell thickness ~ 4 lt-yr, and total mass > 1.2 
Mq, assuming a gas-to-dust ratio ~ 500. Recall that 
SN 1987A is also surrounded by a complex CSM consist- 
ing of three illuminated rings and a dusty peanut-shaped 
structure, located ~ 1 — 2 lt-yr from the progenitor. If SN 
1987A is indicative of the mass-loss structures expected 
around high-mass stars, then a large wind-blown bubble 



(representing main-sequence and/or early red-supergiant 
winds) is expected to be accompanied by inner circum- 
stellar structures that were fo rmed shortly before core 
collapse. Indeed, Dwck (1983]) invoked thermal echoes 
from such inner material to explain early-time IR ex- 
cesses from SN 1980K. 

The deduced structure could also be consistent with an 
evolved version of the type of shell found around Lumi- 
nous Blue Variables (LBVs) such as AG Car, He 3-591 
and Sk —69° 279, with respective shell dimensio ns of 
6.5x 4.5 lt-yr, 1.6x1.6 lt-yr, and 20.2x14.7 lt-yr. (jWeisI 
2011). The nebulae around the massive central stars 
of the AG Car and He 3-591 nebulae have been esti- 
mated to ha ve cool dust masse s of 0.25 M Q and M Q , 
respectively (jVoors et al.l [2000). The mass loss rate 
of AG Car has been estimated to r ange between 1.5 - 
3.7xl0~ 5 Mq yr" 1 (|Barlowi Il99lt IGroh et al.l 120091) . 
comparable to the mass loss rate es timate d for t he im- 
mediate progenitor of SN 1980k by iDwekl (|1983[) . The 
mid-IR spectrum of AG Car is featureless apart from a 
weak 11.3 /im PAH emission band (Voors et al. 2000, Ra- 
jagopal et al. 2007), also consistent with the dust prop- 
erties found here to be needed to match the spectrum of 
the light echo emission around SN 1980k. Regrettably, 
the observations of and constraints from all echoes are 
insufficient to warrant a more thorough morphological 
investigation of these proposed inner and outer circum- 
stellar echoes. 

As discussed in £14.31 an echo will pass out of a spherical 
shell of size t/2 lt-yr in t years, thus it is important to 
note that the best-fit models imply that the echo left 
the CSM shortly after t = 30 yrs. Fig. shows the 
expected V and 4.5 fim light curves for thin (rO = 15, 
6=1) and thick (rO = 13.5, 6 = 5) shells that are 
consistent with the data. Note the precipitous drop in 
flux in the thin shell, as compared to the slow fading 
in thicker shells. Taken on its own, the drop observed 
in the 4.5 /jm datum from t = 30.1 yrs (see also Fig. 
[2]) is consistent with the light curve expected from the 
thin-shell model. It should be emphasized that such a 
drop in the light curve after 30 yrs was not fit during our 
modeling. Obviously, future observations in the optical 
and mid-IR are required to ensure the SN is actually 
fading; the epoch of this fading and its duration will serve 
to disentangle the broad range of acceptable models in 
Fig. El Still, with the data in hand it appears that light 
echoes from a thin shell with inner radius around 14— 
15 lt-yr are an excellent explanation for the multiband 
observations of SN 1980K. 

5. COMPARISON WITH PREVIOUS ANALYSES 

As noted in §1 iFesen fc Beckerl (I1990D . 
Fesen fc Matonickl (|1994D . and IFesen et al.l (|1995l 
19991 ) have reported a slowly-varying complex of 
highly-blueshifted (y r < 7000 km s-1) Ha and [O I] 
AA6300,6363 emission lines. Repres entative spectra 
are re produced at the top of Fig. [8j IFesen fc Beckerl 
(1990) ruled out these high-velocity features as being 
due to a light echo, since such features were not seen 
in early spectra. Instead, they proposed that such 
features might arise from ejecta shocking pre-existing 
circumsteller material. It is therefore important to 
consider whether the observational signatures of such 
an impact are consistent with our proposed light-echo 
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Figure 7. Light curves in V (bottom) and at 4.5 /im (top) of 
echoes off of thin (solid) and thick (dashed) spherical shells that fit 
the observed data, along with the observed data in both wavebands. 

scenario in explaining the last few decades of evolution. 

The only very-young supernova remnant that has been 
visually confirmed to have ejecta-CSM impact is SN 
1987A, in which the forward blast has been lighting up 
"hot spots" around t he in ner equatorial ring or "ER" 
(e.g. iSugerman et al.ll200l iFrance et alJl2010( ). Repre- 
sentative spectra from roughly the same epochs as the 
very-late-time spectra from SN 1980K are shown at the 
bottom of Fig. E The day 2873 spectrum is an HST 
Cosmic-Origins Spectrograph observation of the entire 
inner ring, first hot-spot and ejecta taken on 1995 Jan 
08, and the day 8723 spectrum is of the ring and a few 
bright hot spots taken with the Space Telescope Imag- 
ing Spectrograph on 2009 Oct 21. Both spectra were 
extracted from pipeline-calibrated data taken from the 
HST archive. Note the narrow Ha, [N II] and [O I] lines 
arising from the unshocked ring material. Whether in- 
cluding the entire ring and impact system (day 2873), or 
only a subset (day 8273), these spectra do not show the 
very- high blueshift emission lines seen i n SN 1980K . 

In the plane-parallel approximation ([Sgro| [l975f). the 
transmitted shock velocity vt when impacting media with 
overdensity 5 is 




where v e is the ejecta shock speed and j3 — 6(1 + 
1.701S- 1 / 2 - 0.4018a- 1 + 0.2274<T 3 / 2 - 0.0874(5- 2 )- 2 
ijBorkowski et all 119971) measures the pressure differen- 
tial pre- and post-shock. While the SN 1987A ejecta 
are moving at ~ 3500 km s , the shocks transmitted 
into the ER travel at significantly lower speeds (~ 200 
km s" 1 ), producin g lines that are far narrower and 
less doppler-shifted ([Michael et all2000HPun et al.ll2002D 
than those in SN 1980K. It is therefore unclear whether 
the high- velocity features in SN 1980K can be adequately 
explained by an impact with slow-moving CSM. 

Nonetheless, if this impact is indeed underway, then 
dust present within the CSM will be heated either col- 
lisionally or ra diatively within the sh ocked and/or post- 
shocked gas. iBouchet et all ([20 06) have investigated 
these mechanisms for thermal emission from the hot 
spots in SN 1987A, and find that the dust is most-likely 
heated by radiation from post-shocked gas, has a den- 
sity ~ 10~ 8 cm -3 , and is heated to ~ 150 K. The SN 
1980K CSM models from iDwekl (fl983l ) yield dust and 
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Figure 8. Spectra of SN 1980K (top four), SN 1998S (middle 
two), and SN 1987A (bottom three) at the dates indicated, and 
scaled to arbitrary flux. In grey under the SN 1980K spectrum at 
day 6210 is plotted the estimated fluence of the SN from Fig. [4ji 
as it would appear in an echo when scattered by MRN dust. The 
dashed thin grey lines show the rest wavelengths for Ha and the 
[O I] AA6300,6363 doublet. 

gas densities similar to those i n the ER of SN 19 87A, 
while spectral analyses from e.g. iFesen et all (]1999f ) sug- 
gest the ejecta could be expanding up to ~ 5000 km s _1 . 
Since post-shock gas temperature scales as u 2 (see Eq. 
(fIU[) below), dust that is radiativ ely heated by that gas 
will have temperature Td oc y/vt ([Borkowski et al.lll997t 
IBouchet et aT]|20Ml . We can therefore approximate the 
ejecta-CSM impact as heating dust to < 300 K, which 
is close to the range of warm dust found in our two- 
temperature fits in g3?3l A single "skin" of 300 K MRN 
dust t hat is 10 16 c m thick at the ~ 3 x 10 17 cm radius 
of the IDwekl (|1983D model emits up to an order of mag- 
nitude more flux at 8 /xm than has been observed. I f 
one adopts the two-shock model from lPun et all (|2002f) . 
in which slow and fast shocks are transmitted into the 
CSM based on the geometry and density of the impacted 
medium, and if one makes the reasonable assumption 
that the CSM is not a complete spheroid (i.e. thermal 
emission is not coming from all Att sr) then an ejecta 
impact could explain the average mid-IR SED, provided 
the slow shock component has a velocity < 25% that of 
the fast shocks. 

A complete investigation of the radiation expected 
from ejecta-CSM impact requires time-dependent gas 
and dust radiative-transfer shock modeling, which is out- 
side the scope of this paper. However we can compare 
the heating and cooling timescales to estimate the over- 
all energetic evolution of s uch a syst e m. A gain taking 
SN 1987A as our example, iPun et all (|2002j ) found that 
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a forward blast traveling at ~ 3500 km s -1 from an H II 
region with n ~ 150 cm -3 into CSM with n > 10 4 cm -3 
will produce shocks traveling at v t ~ 250 km s _1 and 
heating post-shock gas to an equilibrium temperature of 
~ 10 6 K. Using their cooling-time approximation 



2 x 10 amu cm 



250 km s~ 



( 10 ) 

where no is the pre-shock density of the impacted ma- 
terial, and [i a is the mean atomic weight, we expect ra- 
diative cooling times of order 1-2 yr. With similar den- 
sities and transmitted shock velocities (preceding para- 
graph) we expect the cooling timescales of shocked CSM 
in SN 1980K to be of this same order. To explain the 
high-velocity lines and thermal emission (i.e. seen 10— 
30 yr after outburst), the ejecta-CSM shocks would have 
characteristic ages t s h oc k 3> t coo i, meaning the observed 
flux levels scale with the surface area of material being 
shocked. For a complete or partial shell, this would in- 
crease with time, which is inconsistent with the observed 
fading across all wavebands. It is certainly to be expected 
that the ejecta from SN 1980K are impacting inner CSM, 
however it appears that this mechanism can account for 
(at most) a small fraction of the observed mid-IR flux, 
and may be predominantly responsible for the presence 
of the narrow systemic emission lines. 

The shape and evolution of Ha emission lines from 
SN 1998S are also indicative of ejecta-CSM interaction 
(jPozzo et alJ2 004): two representative spectra are shown 
in the middle of Fig. [5] Note the early appearance of 
a triple-peaked Ha profile, believ ed to arise from im - 
pact of ejecta with an edge-on ring (IGerardv et alJl2000h . 
By day 658, the red wing of this profile had disap- 
peared, likely due to preferential extinction of emission 
from receding ejecta by dust that had condensed within 
it. This a ffect was also see n in the [O I] lines in SN 
1987A by (|Lucv et all 11989ft . the spectrum of which is 
shown for day 736 in Fig. [3] At later times, the emer- 
gence of a strong, central peak, slowing of the blueshifted 
line, and lack of any red emission can be interpreted by 
"clumpy wind" model of iChugai fc Danzigerl (fl99l in 
which intermediate-width spectral features seen at late 
times arise from shocked clumps of dense wind. 

Although SN 1998S is a different Type (Iln) of SN, The 
Ha profile in its late spectrum bears enough resemblance 
to those of the Ha and [O I] lines in the very-late spectra 
of SN 1980K to suggest that the same mechanism may 
be at play. One concern with invoking ejecta shocking a 
clumpy progenitor wind is that slow-moving clumps must 
be accelerated to high velocities to explain the observa- 
tions, and simulations of ejecta-clump impacts to model 
the hot spots in SN 1987A result in post-shock gas mov- 
ing at hundreds of km s 1 rather than many thousands, 
as needed to match the systemic speeds of the high- 
velocity emission lines. We question whether these high- 
velocity features may instead arise from reverse shocks 
traveling through high- velocity and high-density clumps 
of metal-rich ejecta. The se have long been predic ted to 
form in 2-D simulations ( Chev alier fc Klein IT978). and 
recently observed in fully 3-D simulations of a symmetric 
core-collapse explosions (|Hammer et al.l l2010). Further- 
more, such clumps are common in evolved remnants such 



as Cas A dFesen et al.ll201lD. N63A dWarren et al .1 12003ft, 
and G292.0+1.8 (jGhavamian et al.ll2005ft ~ 

As above, the time-dependent modeling of such a sys- 
tem is outside the scope of this paper, however it is 
instructive to considering whether emission from such 
clumps is expec ted to evolve over the observed multi- 
year timescales (|Fesen et al.l I1999T) . We ad o pt the the 
homologous-expansion model of iChevalierl (|1982ft . in 
which the size R of a homologously expanding medium 
scales as R cx f( n_3 )/( n_s ) where n and s are the density 
power-law indices of the expanding and stationary media, 
respectively. Making the typical assumption that n = 9 
and s = 2 for ejecta expanding into a low-density CSM, 
the size of an ejecta clum p will grow roughly as t & / 7 . In 
the lHammer et afl ()2010[ ) simulation, oxygen-rich clumps 
(which would produce the observed high-velocity [O I] 
lines) have sizes around ~2x 10 11 cm and gas densities 
~ 10 16 cm -3 about 9000 s after core collapse. At an age 
of ~ 300 yrs (roughly that of Cas A) , these clumps would 
grow to R ~ 3 x 10 16 cm and n ~ 100 cm~ 3 , which are 
in excellent agreement with size scales of 1 — 2 x 10 16 cm 
and de nsities of 100 — 300 cm" 3 for metal-rich kno ts in 
Cas A (|Fesen et al.ll2001l : iDocenko fc Sunvaevl[20ToT ). 

At ages of 15 and 30 yr, such clumps would have sizes 
R ~ 1.5 x 10 15 and ~4x 10 15 cm, and densities n ~ 10 4 
and ~ 2 x 10 3 cm -3 , respectively, which correspond 
roughly to the epochs of spectroscopic observation of 
high- velocity lines, and mid-IR photometry, respectively. 
An oxygen-dominated clump (pa ~ 18) traveling out- 
ward at thousands of km s _1 will have radiative cooling 
times C 1 yr using Eq. (|10j) . which is much shorter than 
the crossing times of 1-4 yr for the expected range of 
transmitted shock velocities. As above, a shocked clump 
will emit over its crossing time, which corresponds well 
with the timescales during which the [O I] lines have been 
observed to change. 

An oxygen-rich clump with the properties above will 
have a total mass ~ 10 -5 M Q . If silicate dust within 
this clump forms with an efficiency of 1%, its thermal 
energy output at a few hundred K will be over an order- 
of-magnitude fainter than the observed mid-IR fluxes. 
However this assumes that the entire clump is radiating, 
and given the short cooling times, it is likely that the 
mid-IR emission from a given clump at a given time will 
be significantly less. We cannot rule out the possibil- 
ity that a number of such clumps with a variety of dif- 
ferent pre-shock densities and geometries could provide 
the temperature range needed to fit the thermal SEDs. 
As with the ejecta/CSM impact however, one expects 
that the number of emitting loci will increase with time 
as more clumps are shocked, in contrast to our obser- 
vations of a slow mid-IR decline. We posit that while 
high- velocity clumps may explain the observed emission- 
line structures, they contribute only a small fraction of 
the mid-IR emission. 

Behind the day 6210 spectrum of SN 1980K (Fig. [8j we 
have plotted the fluence of SN 1980K (g2J) as scattered 
by the best-fit thin dust shell. Note that this echo spec- 
trum can account for most of the very broad Ha profile, 
as well as the faint continua present in all the late-time 
spectra. Thus we believe that the proposed scattered 
and thermal echoes, in combination with high- velocity, 
oxygen-rich ejecta clumps, are a reasonable explanation 
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for the observed photometric and spectroscopic very-late- 
time behavior of SN 1980K. 

6. CONCLUSIONS 
Below we summarize the main results of this paper. 

1. We have observed SN 1980K in optical and mid- 
IR wavebands from 23 to 30 years after outburst. 
The SN was not detected in the far-IR bands of 
Herschel. 

2. The optical light curves are consistent with fading 
that has been previously reported, and the mid-IR 
light curves all show the same trend. 

3. The instantaneous dust mass needed to explain the 
mid-IR SEDs is > 10~ 3 Mq, with the presence of 
as much as a few M Q of cold dust possible given 
the large upper limits in the far-IR. 

4. We have presented a new method for computing the 
effective optic al dep th of a scattering and absorbing 
medium (Eq. [A12] ). and an empirical correction 
that is independen t of t he particular grain species 
or properties (Eq. |A14j ^1 . 

5. When used to characterize the line-of-sight extinc- 
tion to a particular location within a homogenous 
dusty medium, these expressions for the effective 
optical depth in the equations of scattering (Eq. [5] ) 
and radiative balance (Eq. [7] ) yield dust tempera- 
tures and emerging scattered and thermal spectra 
that are correct to within ~ 10% for optically-thin, 
and to within a factor of two for optically-thick me- 
dia, when compared to solutions using the Monte- 
Carlo radiative-transfer code MOCASSIN. In prac- 
tice, this "QuickSAND" runs up to a few hundred 
times faster than a traditional Monte-Carlo dust 
radiative-transfer code. 

6. We use the QuickSAND model to compute the 
thermal and scattered echo from a SN over a 
large parameter space of dust composition, density 
power laws, and inner and outer shell radii, and 
find that the observed optical and mid-IR SEDs 
are consistent with echoes from a thin CSM shell 
of inner radius ~ 14 lt-yr and mass < 0.02 Mq 
of carbon-rich dust. This is consistent with the 
size expected for the contact discontinuity between 
a massive star's winds and the ambient ISM, and 
with the size, composition and mass of the contact 
discontinuity found around SN 1987A. 

7. Previous observations of SN 1980K have revealed 
highly-blueshifted [O I] and Ha lines, interpreted as 
arising in impact of the SN blast with pre-existing 
CSM. "Back of the envelope" considerations sug- 
gest that these lines may arise instead in metal rich 
clumps of ejecta, while the ejecta-CSM impact may 
give rise instead to narrower lines observed at their 
rest wavelengths, and that neither mechanism is 
likely to produce sufficient thermal emission to ex- 
plain the observed mid-IR SEDs. 



The quality of the light-echo modeling presented here 
is not exclusive to MRN dust for a spherical shell. For 
example, limiting the largest grains to sizes < 0.1 fixa 
(as suggested by ligh t-echo studies of inner CSM by 
ISugerman et al.ll2005l ) yields acceptable fits with nearly 
the same geometry but a total dust mass of up to 0.06 
Mq. Certainly one could study an enormous range of 
grain-size distributions, optical constants and shell ge- 
ometries (i.e. ellipsoids, disks, hourglasses, etc.) and we 
expect that many of these would yield roughly similar 
fits. Given that this system is unresolved and sparsely 
sampled in the optical, we must accept the ambiguous re- 
sult that the total dust mass in the proposed spheroidal 
wind-blown bubble is < 0.1 Mq, until another decade 
of contemporaneous optical and mid-IR observations are 
available. 

In reconstru ct ing the CSM of SN 1987A, 
ISugerman etaLl ((2005) found a total of < 0.02 
Mq of dust (based on the ambiguity of the gas-to-dust 
mass ratio), with that budget dominated by the large 
contact-discontinuity. While close to the quasi upper- 
limit of < 0.1 Mq of dust in the same shell around SN 
1980K, both values are under the total mass of 0.3 — 3 
Mq of dust in the SN 1987A remnant inferre d by recent 
Herschel observations iMatsuura et al.ll201lD . As shown 
in Tabled a large amount of cool dust can reside in the 
SN 1980K system and remain undetected by Herschel at 
the depth of the 2010 observations. It remains unclear 
whether the total dust mass present in this remnant is 
consistent with that in SN 1987A. Significantly deeper 
Herschel observations are needed of SN 1980K to place 
more stringent limits on this value. 

The early SN 1980K spectra used in this work are 
part of the Asiago-Padova Supernova Archive. We 
thank Annop Wongwathanarat and the authors of 
IHammer et al.l (|2010D for generously providing the data 
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referee, Fernando Patat, for his critical reading and in- 
sightful feedback on the manuscript, as well as bringing 
the ASA spectra of SN 1980K to our attention. Support 
for B.E.K.S. for this work was provided by Hubble Space 
Telescope award GO-10607. This work is based in part 
on data obtained with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with NASA; 
with the NASA/European Space Agency Hubble Space 
Telescope, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. (AURA), under 
NASA contract NAS 5-26555; with the Gemini Obser- 
vatory, which is operated by AURA under a cooperative 
agreement with the NSF on behalf of the Gemini partner- 
ship; and with Herschel, an ESA space observatory with 
science instruments provided by European-led Principal 
Investigator consortia and with important participation 
from NASA. IRAF is distributed by the National Optical 
Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science 
Foundation. 
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APPENDIX 

A. SEMI-ANALYTIC DUST RADIATIVE TRANSFER 
A.l. Isotropic Scattering 

Following the arguments in iRvbickv fc Lightmanl (|1985), consider a homogenous, dusty medium, characterized by 
absorption and scattering cross sections C a bs and C sca , respectively. Traditionally, the optical depth r to some distance 
L is defined as r = L/l where 

Cabs "i" Csca 

is the mean- free path inside that medium. Since (C a b s + C sca ) is the total extinction cross section C ex t, it follows that 
the fractional absorption 

^abs t ^sca 

which also gives the probability that a photon path will end with absorption. If the physical size L of the medium 
is large compared to the mean-free path, then the average number N of mean-free paths that the photon will travel 
before being absorbed is 

N=l/e (A3) 

since e is the probability of a given step ending in absorption. 

The path of a photon in that medium can be expressed as the vector sum R = $D i=1 ^ of the N scatters taken 
before absorption (or escape). For isotropic scattering, a photon's path is a random walk, with the well-known results 

that the average final position (R) of many photons will be zero, while the average displacement of the photons from 



their starting point is the RMS distance I* — y (R 2 ), which expands out to 

N N N 

n=Y,tf)+2Yl E ( A4 ) 

2—1 2—1 

Since the scattering is assumed to be isotropic, the first term (rf) ps I 2 while the cosines in the cross terms average to 
zero, yielding the expected result 2* = ly/N. Combining this with Eqs. (|A1 |> — (|A3 JI yields I 2 — l 2 /e or 

I* = [C a bs(C a bs + Csca,)] ^ ■ (A5) 

The optical thickness of a medium is approximately the number of mean- free paths required to escape, i.e. r = L/l 
as defined above. For an optically thick medium where the probability of absorption is non trivial, one must redefine 
optical thickness as the ratio of the physical size L of the medium to the average distance 1* between creation and 
destruction of a photon. As such, 1* is often called the "effective mean path" and the resulting effective optical depth 

r* = - (A6) 

(A7) 



T 



where the optical depth due to absorption, scattering, and extinction are r a b s = LC a b s , t sc = LC sca , and r = r a b s 
respectively. Defining the scattering albedo u> = t sc /t, a little algebra yields 



r* = tVI - oj (A8) 

which is equivalent to Eq. (14) of VD99. 

A. 2. Non-Isotropic Scattering 

If one relaxes the assumption of isotropic scattering, a more generalized expression can be derived by noting that the 
scattering anisotropy parameter g = (cos 6), where 9 is the scattering angle between incoming and outgoing photon 
traje ctories. For isotropic scattering, g — 0, while g — > 1 for pure forward scattering. If g ^ 0, the cross terms in Eq. 
(|A4p . (r*j- • ~ gl 2 , and by recursion, (r*j ■ fj) w gi~ % l 2 . Since g < 1, the cross terms form a convergent geometric 

sequence such that 

i=i j=i+i ^ y ' 

Inserting this back into Eq. (|A4I) yields 

T = (MO) 
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Figure 9. Selected results from testin g of th e effective optical depth, showing r* — r e g for different values of the para meter triplet ui, g and 
r. The top row computes t* with Eq. I)A12|| . The bottom row includes the empirical correction term St* (Eq. |A14| ). Positive (negative) 
contours are plotted with solid (dashed) lines at the levels indicated. 
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Figure 10. Comparison of our semi-analytic dust radiative-transfer model (QuickSAND) to MOCASSIN for the SN 1980K UV (top row) 
and optical (bottom row) fluences illuminating spherical shells of optical depth listed at the top left corner of each panel. The SED from 
MOCASSIN is shown in thick grey, and QuickSAND in thin black in the top half of each panel pair. The bottom half of each pair shows 
the ratio of the two models. 



where 

N - g (2 + gN - 2g N ) 
N.(g) = \ g _ l? ' ( A11 ) 

and N is defined as in Eq. (|A"3]l . Note that the limit N*(g -> 1) = N 2 . Finally, substituting Eq. (|XT0)l into Eq. |MJ 
yields the more general expression for effective optical depth 

rM 'vm- (A12) 

Note that Eq. (|A12[) reduces to Eq. (IA7|) in for isotropic (g = 0) scattering, while for fully forward-scattering particles 
(g = 1), T»(g) « T /N which is equivalent to Eq. (15) of VD99. 
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A. 3. Testing 

We consider the dust radiative-transfer code MOCASSIN (jErcolano et al.l l2005) to be the benchmark against which 
the effective optical depth r* can be tested. MOCASSIN self-consistently computes the radiative-transfer of a given 
energy spectrum Fq(X) through a dusty medium composed of user-specified grain properties within a Cartesian grid 
of cells by following the paths of energy packets (i.e. photons) until the emergent flux has converged to a specified 
tolerance. Since photons follow random paths, the emergent spectrum F(X) has an inherent variance that depends on 
the number of photons and the convergence tolerance. 

By comparing the emergent and input (or unattenuated) fluxes, the effective optical depth of a given model is defined 

as 

« = lnf§. (A13) 

A grid of models was constructed with ui and g ranging from 0. to 1. and optical depth r from 0.1 to 15, allowing us to 
compare t» to r e g for each triplet of parameters. Rather than using actual dust-scattering properties, each model used 
uniform values of C a b s , C'sca and g for a variety of wavelengths within a spherical shell of inner and outer radii 10 
and 10 15 cm to match the parameter triplet (u>, r, g) being tested. Using a uniform input spectrum, r e ff is computed 
as the average across all wavelengths, providing a measurement that is more robust to the inherent variance of the 
Monte-Carlo process. Selected results for r = 1, 5 and 10 are shown in Fig. [9] The quality of the approximation t* 
are roughly consistent with those presented in VD99. The most substantial difference is that the residuals r* — T e g are 
well-described by an elliptical gaussian, which is approximately modeled by the empirical expression 

, _ / (0.219r- 0.687)e^ 41n2 ("/°- 75 ) 2 if r < 5.25 , 

* 1 7.5 X l -3 r 2.5 e -41„2(p/0.75) 2 &]se ^> 

where p = (^{g - g ) 2 + {uj - ui ) 2 y 1/2 , g = 0.3 (l + e" 7 / 3 ) and uj = 0.585 (l + e^'*). The bottom row of Fig. [9] 
shows the result of adding this correction term to t» , yielding a very accurate approximation of T e g over a very broad 
range of dust properties in the optically thick regime. 

A.4. The QuickSAND Model 

Replacing the optical depth r in Eqs. (UJ), ([7]), and ((HJ) with r* 4- St* from Eqs. (|Al2j) and (|A14|) . we can compute 
the dust temperature and emergent spectra from every point within a homogenous (i.e. non-clumpy) medium. We call 
this the "QuickSAND" (Quick Semi-ANalytic Dust) radiative-transfer model. 

We benchmark tested the QuickSAND model against MOCASSIN by illuminating a spherical shell of inner and outer 
radii of 10 and 15 ly, respectively, composed of an even mixture of carbonaceous and silicate MRN dust ( §3.3p with 
constant density and illuminated by the UV and optical fluences in Fig. |4] The MOCASSIN models use a reasonably 
dense cubic grid of 61 cells in each dimension and 128 million energy packets, both of which are needed to ensure that 
sufficient photons are scattered into the line of sight to provide adequate signal-to-noise in optical scattering and the 
far-IR. The QuickSAND models sample the shell at roughly double the spatial resolution as MOCASSIN. 

Resulting benchmark SEDs from both models are shown in Fig. [10] for Ty = 0.01, 0.1, 1., and 10. The QuickSAND 
model provides an excellent approximation to the MOCASSIN model at low and intermediate optical depths, and is 
correct to within a factor of two in optically-thick media. For comparison purposes, the MOCASSIN model running 
in parallel on eight processors requires roughly 30 times longer to converge than the QuickSAND model running on a 
single processor. 
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